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Abstract 

We discuss the Poisson structures, Lax matrices, r-matrices, bi-hamiltonian structures, 
the variables of separation and other attributes of the modern theory of dynamical systems 
in application to the integrable Euler top and to the nonholonomic Chaplygin ball. 



1 Introduction 

The main aim of this paper is to prove that the integrable Euler top and the nonholonomic 
Chaplygin ball are very similar dynamical systems like birds of a feather flock together. Thus, 
on example of these twins, we want to show how all the machinery developed for integrable 
systems can be carried to the theory of solvable nonholonomic systems. 

The integrable Euler case of rigid body motion with the fixed center of mass (the Euler top) 
is relatively simple in the sense that its equations of motion do not linearize on Abelian surfaces, 
but on the elliptic curves. Of course, this does not make the Euler top entirely trivial [TB]- A 
classical description of the Euler top can be found in any textbook on classical mechanics, see, 
for instance, [HE1E]. 

The nonholonomic Chaplygin ball is that of a dynamically balanced 3-dimensional 
ball that rolls on a horizontal table without slipping or sliding. 'Dynamically balanced' means 
that the geometric center coincides with the center of mass but the mass distribution is not 
assumed to be homogeneous. Because of the roughness of the table this ball cannot slip, but 
it can turn about the vertical axis without violating the constraints. There is a large body 
of literature dedicated to the Chaplygin ball, including the study of its generalizations. Sec 
[H [HI 121 H31 HH [IHl I2SJ HSl ISI]- Of course, this list, as well as the bibliography of the present 
paper, is by far incomplete. 

Section 2 starts by collecting some definitions and facts about Euler top and Chaplygin 
ball. In Section 3 we will consider Poisson structures associated with the Chaplygin ball as de- 
formations of the similar standard Poisson structures for the Euler top. Section 4 contains our 
main results about separability at zero level of the cyclic integral of motion. In bi-hamiltonian 
geometry separability is an invariant geometric property of the distribution defined by mutually 
commuting independent integrals of motion. In fact, there is neither Hamilton- Jacobi equa- 
tion, nor time which describes only partial parametrization of geometric objects. We want to 
show how those standard bi-Hamiltonian geometric methods may be directly applied to the 
nonholonomic Chaplygin ball. 

In the second part of the paper we will discuss various deformations of the well-known 
integrable Hamiltonian systems which may be treated as generalizations of the Chaplygin ball. 
In Section 5 we will briefly consider deformations of integrable systems on cotangent bundles 
of the Ricmannian manifolds and underline that the main problem is the change of time which 
transfers a purely mathematical construction to the sensible physical model. For all these 
systems the deformations of initial Poisson bracket is trivial. Finally in Section 6 we will give 
some examples of similar deformations on Lie algebra e(3), when the deformation of the initial 
Poisson bracket is nontrivial. 
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2 Equations of motion 

Let 7 = (71,72,73) and M = (Mi, M2, M3) be the two vectors of coordinates and momenta, 
respectively. We postulate that they satisfy to the following differential equations 

M = M x uj , 7 = 7 x uj . (1) 

For any vector function u) on the dynamical variables x = 7, M these equations in A4 = I 3 x I 3 
have the following integrals of motion 

Hi = (7,7), n 2 = (%M), H 3 = (M,M). (2) 

Six differential equations can be solved in quadratures if we know four integrals and the Jacobi 
multiplier |18j . So, we want to add some additional integral to the known integrals T-L\, I-L2 and 
H3 ([!]) and calculate the desired multiplier. 

If we assume the existence of the following additional integral of motion 

Hi = (M,uj) (3) 

one gets 

^ = [M x to, uj) + (A/, to) = (M, w) = , 

it means that the derivative Cj has to be perpendicular to M. Below we stint ourselves by 
integrals © with 

u = A X M . (4) 

In generic A x is a matrix depending on variables x — (7, M), which has to satisfy to the 
equation 

(M,w) = (AjM, M x A X M) + (M,A X M) = . 
This equation can be replaced by the particular system of very simple equations 

Aj=A x , and (M,A x M)=0, (5) 

which has constant solution associated with the Euler top 




A x = A, A= I * / j () j . a, I! «i) 

and solution associated with the Chaplygin ball 

A x = A d , where A d = A + dg(-y) A7 ® 7 A , (7) 

A d = g(7) 2 A (7 ® /3 + /3 ® 7) A, /3 = (7- d( 7 ,7)A 7 ) x AM. (8) 
Here A is given by (j6j) and function g(7) is equal to 

so that Ad goes to A at d —> 0. 

In the Euler-Poisson case, M is the vector of the kinetic momentum, uj is the angular 
velocity and 7 is the unit Poisson vector [3J [7] . All these vectors are expressed in the so-called 
body frame, its axes coincide with the principal inertia axes so that the corresponding tensor 
of inertia reads as 

J = A 1 = ( J 2 ] , J k <J t + J 3 . (10) 
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In the Chaplygin case M is the vector of angular momentum of the ball with respect to 
the contact point, u> is the angular velocity vector of the rolling ball and 7 is the unit normal 
vector to the plane at the contact point. As above all these vectors are expressed in the body 
frame firmly attached to the ball. If the mass, the inertia tensor and radius of the rolling ball 
are denoted by to, J and a, then matrix A in ([7]) is equal to 



A = 



V 



1 



J 1 + d 








J 2 + d 




d = 



(11) 



■h + d / 



At d — > equations of motion ([T]) and integrals of motion T~Lk for the Chaplygin ball coincide 
with equations and integrals for the Euler top. Of course, it is purely mathematical and non- 
physical limit because Ji — ¥ as d — > 0. So, from mathematical point of view, we can say that 
Chaplygin ball is a deformation of the Euler top with respect to parameter d. 

In the Euler case matrix A and integral of motion 7^4 are different from the matrix A and 
integral of motion in the Chaplygin case. Nevertheless, for the brevity, we will use common 
notations A and in both cases where it will not cause any confusion. 

Thus, we obtain the fourth integral of motion W4 ((4]) for the six equations (JTJ at w = 
AM and u> — A^M. We proceed by showing that these dynamical systems are solvable in 
quadratures in framework of the Euler- Jacobi last multiplier theory [18j . By definition, the 
Jacobi multiplier fi(x) of (fTJ) is a function on dynamical variables x = 7, M , which has to 
satisfy to the equation 



6 d 3 / d d \ 



= 0. 



For the solution A © this equation is trivial 

ii(x) — and \i(x) — /1 = c, 



c € 



but for Arf one gets 

2g(7)AO) - K x Mi) = o : 



and 



(12) 



(13) 



According to [TH] the Jacobi's multiplier is some nontrivial function in the case of constrained 
systems only. The integrability conditions of the nonholonomic systems formulated by Kozlov 
[26] include the preservation of measure related with the Jacobi multiplier. 

There are many other solutions of the system ([5]), see review |S]. For instance, solution 

A / = fil) 7 ® 7 . 

depending on arbitrary function f(j) is associated with multiplier fi(x) = 1. Solutions associ- 
ated with nontrivial Jacobi multiplier are given by linear in variables 7 matrices 



A abc = A + B( 7 ® c + c <g> 7)B~ 



(14) 



which satisfy ([5]), if we impose various restrictions on the numerical entries of matrices A.B 
and vector c. 

To sum up, we can easily get a lot of additional integrals (U) and the corresponding 
Jacobi multipliers of the equations ([1} and, therefore, we can solve these differential equations 
in quadratures without any notion of the Hamilton structure, integrability by Liouville, the 
Poisson structure, the Lax matrices, classical r-matrices etc. 

However, this additional and in some sense redundant information can be useful in various 
applications, such as the perturbation theory, the quantization theory and so on. Below we 
reconstruct this information starting with only integrals. 
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3 The Poisson brackets. 



In this section our aim is to calculate the Poisson brackets for the given models without any 
assumptions on underlying Hamiltonian or conformally Hamiltonian structures of the equations 
of motion [6l [8] . We will calculate the desired Poisson brackets only assuming that the foliation 
Hi = a.i is a direct sum of symplectic and lagrangian foliations. 

Let us consider the manifold M, dim M. — n, with coordinates x — (xi, . . . ,x n ). The 
Jacobi last multiplier theorem [18] ensures that n equations 

dx ■ 

— = Xi(xi,...,x n ) , i = l,...,n. (15) 

are solvable in quadratures if we have n — 2 functionally independent integrals of motion 
and the Jacobi multiplier [i. 

Let us suppose that M. be a Poisson manifold endowed with a Poisson bivector P, so that 
the corresponding Poisson bracket reads as 

In our examples, all the symplectic foliation associated with the Poisson bivector P is rather 
regular. Moreover, all the leaves are affine hyperplanes of codimension k, which are the level 
sets of k globally defined independent Casimir functions Cj 

PdCj=0, j = l,...,k, 

If the symplectic foliation associated with P is rather regular and n = 2m + k, then we can 
determine the set of Hamiltonian systems on Ai 

^- = {Hj,Xi}, i = l,...,n, j = l,...,n-2. (16) 

The invariant volume form for all these Hamiltonian flows may be formally expressed as a 
square of symplectic form on symplectic leaves 

v = ft m , where ft = P" 1 , 

and P c is restriction of P on symplectic leaves. 

Furthermore if one of these Hamiltonian systems (|16p is integrable by Liouville we can 
compare the two integrable systems (TT5j) and (TT6)) using change of time 

t->t s 

and try to extract useful information about the former system from the known properties of 
the latter Hamiltonian system. 

According to the Liouville theorem, for the given Hamiltonian T-Lj equations of motion (|16[) 
are integrable in quadratures on the symplectic leaves, if we have m functionally independent 
integrals of motion Hi = %j , H2 ■ ■ ■ , H m in the involution 

{Hi,Hj} = 0, i,j =l,...,m. 

We can identify all the integrals of the equations (|T5l) with all the integrals of ([16]) 

(Hi, ■ ■ ■ Tin— 2) ~ (Hi, . . . , H m ; Ci, . . . , Cfc) 

only at m = 2, because n = 2m + k and n — 2 — m + k. 

Regularity of symplectic foliation is closely related to the existence of the Jacobi multiplier 
[i. One of the global invariants in Poisson geometry is a modular class. It is an obstruction 
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to the existence of measure in Ai which is invariant under all hamiltonian flows [24j |47l [48] . 
For the manifold A4 endowed with a Poisson bivector P, its modular class is an element of the 
first Poisson cohomology group. In Section 3 we discuss some elements of the second Poisson 
cohomology group and the corresponding Poisson bivectors P' compatible with P, which allows 
us to get variables of separation without any additional information. 

In generic case we can identify only one integral Hi = Hk and consider not only two- 
dimensional systems, but other systems as well. In this case the integrals Hj could be generators 
of some algebra of integrals with respect to the bracket {.,.}, see theory of superintegrable or 
noncommutative integrable systems [16] . 

In our case integrable by Euler-Jacobi equations of motion (flj with integrals Hi, . . . ,Ha 
and multiplier fi are integrable by Liouville after an appropriate change of time, if there is the 
Poisson bivector P such as 

• [P, P] = 0, the Jacobi identity, 



• PdHi = PdHj = 0, only two Casimir functions (17) 

• {Hi,H m } = 0, the involution of the integrals . 

Here [., .] is the Schouten bracket, Hi, ■ ■ . , Hi are the four integrals ([2][3]), and (i,j, I, m) is the 
arbitrary permutation of (1, 2, 3, 4). 

The first equation in (fl7j) guaranties that P is a Poisson bivector. In the second equation 
we define two Casimir elements Hi and Hj of P and assume that rankP = 4. It is a necessary 
condition because by fixing its values one gets the four dimensional symplectic phase space of 
our dynamical system. The third equation provides that the two remaining integrals Hi and 
H m are in involution with respect to the Poisson bracket associated with P. 

If we consider Chaplygin ball as a deformation of the Euler top, it is natural to fix for the 
both systems the same Casimir functions 

i = l,j = 2, Hi=d, H 2 = C 2 (18) 

and integrals of motion 

Z = 3, m = 4, H A = H l , H 3 =H 2 . (19) 

Solutions P associated with another choice of Casimir functions may be obtained from this 
solution by using classical r-matrix theory, see Section 3.3. 



3.1 The linear in momenta Poisson bivectors. 

The system (TIT)) has infinitely many solutions and, therefore, we have to narrow the search 
space and try to get some particular solutions only. In this Section we assume that the entries 
of P are the linear functions in momenta M. 

Proposition 1 In the hypotheses mentioned above the system of equations |J7[ ) has the follow- 
ing linear in momenta solutions: 

Euler case: P = f £ ^ J , (20) 

Chaplygin case: P d = -L= ( ^ J. ] - d^f^f) (M, A 7 ) " " 



Vsh) 



Here g(j) is given by (0) and 






73 


-72 


-73 





7i 


72 


-7i 






r m ) V5UA 1 " v r 



M = 
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The proof consists in the substitution of the linear in momenta anzats for entries of the Poisson 
bivector 

k 

into ([!?)) and in the solution of the resulting algebro-differential equations with respect to 
unknown coefficients fijk{"f)-0 

The Poisson brackets between variables x = 7, M look like 

{Mi,Mj} = EijkMk, {Mi,7j} = eijklk {7i»7j} = 0, (21) 



and 




{M u Mj} d = e ijk ( ~ dy/^fj(M, A 7 ) 7fc ] , 



{Mi,jj} d 



(22) 



Vsh) 



Here e^k is a totally skew-symmetric tensor. 

The first bracket {.,.} is the well studied Lie-Poisson bracket on the Lie algebra e*(3), 
whereas second bracket may be considered as its deformation with respect to parameter d, which 
preserves regular symplectic foliation. The second Poisson bracket (|22|) has been obtained in 
®. 

Formally, in the both cases the Poisson bracket (121II22[) allows us to rewrite the initial 
equations of motion (fTJ) in the Hamiltonian form 

dx dx , rr , 1 „ , .„„, 

^ = JKxjdt^ {H ^ h H= 2 n ^ ^ 

after changing the time variable including the corresponding Jacobi multiplier 

dt' = n(x) dt . (24) 

where fJ.(x) = 1 or fi(x) = /i^, respectively. In the Chaplygin case this transformation has been 
introduced in the Chaplygin work [IT] in order to get the solutions as the functions of the time 
variable. 

It is easy to see that for the Chaplygin ball we can not directly identify this transforma- 
tion with canonical transformations of the extended phase space, which change time and the 
Hamilton function simultaneously [301 IM1 12H I3H1 132] > see also discussion in [5J [T3J I33J • 

Of course, we can get similar Poisson bivectors for other solutions of (O as well. For 
instance, if u) = A a i,M, where 




a + 673 

is a matrix from the family (|T4l . then solution of (IT71) looks like 

T \ aM- 



Pab = 1 



\J ax 3 + b 



r M ) 2(ax 3 + b) 



-(° 




\ 


?)] 



(25) 



In this case the equations of motion (fTJ are integrable by Euler- Jacobi theorem and by the 
Liouville theorem after the corresponding change of time. 
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3.2 Properties of the linear Poisson bivectors. 

The Poisson-Lichnerowicz cohomology of the Poisson manifold was defined in [53] , and it pro- 
vides a good framework to express the deformation and the quantization obstructions, see 

[21 Hang. 

Let us remind some necessary facts from the Poisson geometry. The Poisson manifold A4 
is a smooth (or complex manifold) endowed with the Poisson bivector P fulfilling the Jacobi 
condition 

[P,P]=0 

with respect to the Schouten bracket on the algebra of the multivector fields on AL Other 
Poisson bivector P' is compatible with P if any of its linear combination P + XP' is the Poisson 
bivector, i.e. if 

[P,P']=0. 

Bivectors P' are the 2-cocycles in the Poisson-Lichnerowicz cohomology defined by Poisson 
bivector P on the Poisson manifold M. . They must be compared with the bivectors 

pW=C x (P) => [P,P m ]=0 (26) 

which are the Lie derivative of P along any vector field X on A4. Bivectors p( x ) are 2- 
coboundaries and 2-cocycles simultaneously. However not all cocycles are coboundaries. If X 
is such vector field that the Jacobi condition 

[pW pW] = Q 

is satisfied, then p( x ) (|2"6"|) is called the trivial deformation of the Poisson bivector P. 

Now let us go back to our physical models. The first bivector P (f20|) is the well studied 
Lie-Poisson bivector on the Lie algebra e*(3) of Lie group E(3) of Euclidean motions of R 3 . 
The second bivector Pj (|20p can be treated as its "nonholonomic" deformation related with the 
Chaplygin ball. 



Proposition 2 Bivector Pj H20\) is 2-cocycle in the Poisson-Lichnerowicz cohomology defined 
by canonical Poisson bivector P on e*(3). 

The proof is a straightforward calculation of the Schouten bracket 

[PPi]=0. 

Recall that the Schouten bracket {R, QJ of two bivectors R and Q is trivector and its entries in 
local coordinates x look like 

rn—l ^ ' 

□ 

Proposition 3 In the generic case the Poisson bivector 

P d = ^L=( Q T ^-d^r)(M,Aj)(° Q °) (28) 



Vsh) 



is a nontrivial deformation of the standard Lie-Poisson bivector P , which is a sum of two Lie 
derivatives of P 

P d = C Y (P) + d ^\^ M) L z {P). (29) 
Here entries of the vector fields Y = Y 3 dj and Z = Y] Z J dj are given by 

Y l = Z' l = 0, Y t+3 = j=^=, Z l+3 = ((trA-Id- A) 7) , i= 1,2,3. (30) 

Vg(7) v >l 
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In order to prove first part of this proposition we have to try to solve the following equation 



P d = £x(P), 



(31) 



with respect to unknown X. Recall that in local coordinates the Lie derivative of a bivector P 
along a vector field X reads 



dimM 

E 

fc=i 



X 



k dPV pki dX* 
dx k dx k 



- P 



dXi \ 
dx k J 



Using the modern software for symbolic calculations it is easy to prove that entries of (|31|) form 
the inconsistent system of differential equations. It means that (|31j) is infeasible equation and, 
therefore, cocycle Pd is no coboundary. The second part of the proposition is verified by direct 
calculations. □ 

It is easy to see, that if (7, M) = then Pd = Cy(P) is a trivial deformation with all the 
pleasant mathematical and physical consequences, see [291 H7] and [51 [TT] respectively. 

In the finite-dimensional case local Poisson geometry begins with the splitting theorem, 
which says that in the neighborhood of any point in the Poisson manifold .M, there are coor- 
dinates (qx,. ..,q m ,px,... ,p m , Ci, . . . ,C k ) such as 



P 



V — A — 

% dpi 



\ E ^(°) 



d 



d 



dd A dd 



and ipij (0) = 



So, if the compatible bivectors P and P' have a common set of Casimirs C\, . . . , CV-, we can 
identify the Darboux coordinates (q,p) of P with the Darboux coordinates (q',p r ) of P' and 
obtain the local map cf> : M. A4, which pulls back P' to P. In Section 4.3 we prove that in 
our case this local map Pd — > P can be extended to the global one at C2 = 0. 

If we come back to the general theory, the second Poisson-Lichnerowicz cohomology group 
Ti. 2 p on A4 is precisely the set of bivectors P 1 solving [P, P'J = modulo the solutions of the form 
p( x ) = Cx(P)- We can interpret T-L 2 P as the space of infinitesimal deformations of the Poisson 
structure modulo trivial deformations. We should keep in mind that cohomology reflects the 
topology of the leaf space and the variation in the symplectic structure as one passes from one 
leaf to another [HI Ell Wj\ HE] . 

3.3 The r-matrices 

It is known that equations ([T]) can be rewritten in the Lax form 



dL 

~dt 



[l, n], 



M 



r 

A' 



A G 



(32) 



if we identify (K 3 , x) and (so(3), [., .]) by using a well known isomorphism 







Z3 


z = (z 1 ,z 2l z 3 ) -> Z = 









V Z 2 


~z x 



(33) 



where x is a cross product in IR 3 and [., .] is a matrix commutator in so(3). Another possibility 
is to use 4x4 antisymmetric matrices 



M 






M 3 


-M 2 




-M 3 





M, 





M 2 


-Mi 

















0/ 



ft 






w 3 


-CJ 2 


M 


-CJ 3 





CJl 





UJ-2 


-CJl 

















/ 
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and symmetric matrix 







7i 







72 







73 



71 

72 

73 




In the both cases bilinear Lax matrix L(A) ([3"2"]) belongs to a huge family of Lax matrices 
described in the book [34] , see also [35]. Below we will consider only 3x3 Lax matrices. 

The Lax equation implies that the spectral invariants of the Lax matrix L(A) are conserved 
quantities under the Hamiltonian evolution, but their involutivity and functionally indepen- 
dence must be checked case by case. In a noteworthy paper [3J, Babelon and Viallet showed 
that if all the spectral invariants of the mxm matrix L(A) are in involution with respect to some 
Poisson bracket {., .} on a given phase space, then there is a matrix ri2(A, /i) of order m 2 x m 2 
such that the Poisson brackets between the entries of L are represented in the commutator form 



{L(A),L(^)} fc 



[n 2 (A, n) , L] - [r 2 i (A, /x) , L(/x) ; 



(34) 



1 2 

HereL(A) = L(A) ® Id , L(/i) = Id ® L(/x) and 7*12 (A, fi) is a classical r-matrix and 

»*2i(A,m) = Pri 2 (/i,A) , 

where P is a permutation operator: ILc ® y = y <g> a;, Vx, y £ C m [34] . 

In our case for the Euler top and for the Chaplygin ball we have one Lax matrix L(A) 
and two different Poisson brackets. First bracket {. , .} (|2ip associated with the bivector P ([2"0"]) 
yields the standard r-matrix 



ri2(A,/x) = 



fj, — A 



i=l 



Here S; form a basis in the space of 3 x 3 antisymmetric matrices 




'o -r 

S 2 = [ 

.10 



1 0^ 

S, = 1-1 

0, 



(35) 



Si = 



For the nonholonomic bracket {., .}d ([22]) associated with the bivector P3 ([20]) the r-matrix will 
be a more complicated dynamical r-matrix. 



Proposition 4 The Lax matrix for nonholonomic Chaplygin ball t32\) satisfies the linear r- 
matrix algebra |g^[ ) with the following r-matrix 



ri 2 (A,/x) 



H - 



A 



dA A /iC7)(M J A 7 ) ^Si®Si 



i=l 



(36) 



The proof is straightforward verification of (|34]> for the given Lax matrix.. □ Usually, the 
bilinear Lax matrices (1321) are not very useful in integration of equations of motion. However, 
they are very effective in various geometric applications, for instance see one of the latest 
applications in discrete differential geometry [27] . 

In bi-Hamiltonian geometry we can use bilinear Lax matrices in order to get solutions of 
([T7]) associated with another choice of the Casimir elements. Namely, if ri2 is classical r-matrix 
and tp is intertwining operator, then ri2 o ip is also a classical r-matrix [34] . For a given matrix 
7*12 the r-matrices r\i o (p form a linear Lie pencil, which generates a family of compatible 
Lie-Poisson brackets. For instance, if we take trivial intertwining operators from [34] 



and 



<P2 = H 
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and r- matrix (|35j) . then one gets the following well-known Poisson bivectors for the Euler top 



Pl = ( o o ) ' and P2 = ( *o ) ' 

see [35]. Construction of other brackets on e*(3) associated with classical r- matrices (|35M36I) 
and more sophisticated relations between various classical r-matrices |42) will be discussed in 
forthcoming publication. 



4 Separation of variables at (7, M) = 0. 

Now we address the problem of separation of variables within the theoretical scheme of bi- 
hamiltonian geometry PUSHES]- According to [TT] we can start with the case (7, M) — and 
then reduce the generic case to this particular one. 

In geometry, instead of an additive separation of variables in the partial differential equa- 
tion called the Hamilton- Jacobi equation, we have some invariant geometric property of the 
Lagrangian distribution defined by m independent functions Hi, ... , H m 15, 44l 145) . 

Namely, an m-tuple H\ , . . . , H m of functionally independent functions defines a separable 
foliation on symplectic leaves of M, if there are variables of separation [qi, . . . , q m ,pi, ■ ■ ■ ,p m ) 
and m separated relations of the form 



$i = , with dct 



7^0. (37) 



It simple means that the common level surfaces of H\ , . . . , H rn form rather regular foliation 
of symplectic leaves and every leaf of this lagrangian foliation may be represented as a direct 
product of one-dimensional geometric objects defined by separated relations (|37[) . 

Now let us remind how to get variables of separation in framework of the bi-Hamiltonian 
geometry. The bi-Hamiltonian manifold Ai is a smooth (or complex) manifold endowed with 
two compatible Poisson bivectors P and P' . Dynamical systems on M. with the integrals of 
motion in involution with respect to the both brackets 

{H i ,H j } = {H i ,H j y = 0, i,j = l,...,m, (38) 

are called bi-integrable systems (44] [45] . The bi-involutivity of the integrals of motion (|38|) is 
equivalent to the existence of the control matrix F defined by 

m 

P'dHt = P^FijdHj, i = 1, . . . , m. (39) 
j'=i 

The eigenvalues {q%, . . . , q m ) of F are the coordinates of separation, whereas the suitable nor- 
malized left eigenvectors of F form the generalized Stackel matrix S 

F = S^ 1 diag (qi,..., q m ) S 

which defines the separation relations 

n 

^i=^2S ij (q i ,p i )H j + Ui(qi,Pi) = , i=l,...,m. (40) 
j'=i 

Here the entries of Stackel matrix Sij and the Stackel potentials U, depend only on one pair 
(qi,Pi) of the canonical variables of separation, Casimir functions Cj and, in generic case, on 
the integrals of motion Qjj] [44] [45] . 

In our case the Stackel matrix and the Stackel potentials depend only on variables of 
separation, and it allows us to calculate the canonical transformation from the initial variables 
7, M to the variables of separation explicitly. 
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4.1 Darboux-Nijenhuis coordinates 

In order to get variables of separation according to the general usage of bi-hamiltonian geometry 
firstly we have to calculate the bi-hamiltonian structure for the given systems with integrals 
of motion Hi , H2 (TT51) on manifold A4 with the canonical Poisson bivector P (j2"01) and its 
deformation P&. 

Proposition 5 Let us introduce two vector fields X — ^2X J dj and Xd = ^2XZdj, with the 
following entries: 



X 1 = 0, X l + 3 = 



x d = o, x d + 3 = 



7 x A (7 x M) , i = l,2,3. 

(41) 

7 x A d ( 7 x M) . 



The Poisson bivectors 



P' = C X P and P' d = C Xd P d (42) 

are compatible with the bivectors P and Pd h2U\) respectively. Bivectors J^i?| ) have common 
symplectic leaves 

P'dCi = P' d dCi = 0, P'dC 2 = P' d dC 2 = 0, (43) 
whereas integrals of motion Hi i2 H9\) are in the bi-involution 

{H 1 ,H 2 } = {H 1 ,H 2 }' = 0, {H 1 ,H 2 } d = {H u H 2 }' d = Q (44) 

with respect to the corresponding Poisson brackets at (7, M) = only . 

The proof is a straightforward verification of the corresponding Schouten and Poisson 
brackets in local coordinates. □ Thus, we proved that the Euler top and the nonholonomic 
Chaplygin ball are bi-integrable systems at (7, M) — 0. At second step we have to calculate 
the corresponding control matrices F and Fd defined by For the Euler top we have 

_ / (A v 7, 7 ) \ 

F - \ -(7,7) ((trA.Id-A) 7 , 7 ) ) {A5 > 

and similar to the Chaplygin ball 

_ / (A^ 7 ,7) \ 

Fd -[-(7,7) ((trA d .Id-A rf )7, 7 ) )' m 

Here A v = (det A) A -1 is adjoint or cofactor matrix. 

The Darboux-Nijenhuis coordinates associated with the bivectors (|42|) and control matrices 
(|45H46|) are the roots of their characteristic polynomials 



r(A) = A 2 - (jtrA • Id - A) 7 , 7J A + (7, 7)(A V 7 , 7) = , 

(47) 

r d (X) - A 2 - ((trA d • Id - A d ) 7 , 7) A + (7, 7 )(A^7, 7) - . 

By definition [15] . the Darboux-Nijenhuis variables are canonical with respect to the symplectic 
form f2 associated with first bivector P and put the recursion operator TV = P'P^ 1 in diagonal 
form on symplectic leaves of M. 

Below we prove that these Darboux-Nijenhuis coordinates (|4T1) are the variables of separa- 
tion for the bi-lagrangian foliation defined by integrals Hi 2 (fl9f on symplectic leaves of e*(3) 
fixed by d = 1 and C 2 = 0. 

It is easy to see that at (7, M) = the passage from the Euler top to the nonholonomic 
Chaplygin ball consists of the replacement of the constant matrix A © on dynamical one A^ 
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([7|) in the equations of motion ([T]), in the Hamiltonian Hi = (M, AM) and equations (|41|45|47[) 
only. Similar to geometric quantization theory, simplicity of this deformation is a sequence of 
the equation (|2"5)) 

P d = C Y {P), 

properties of the Lie derivative C and of the vector fields Y (150)) and X4 (|4"Tj) . 



4.2 Elliptic coordinates 

At C 2 = and C\ = 1 we can identify the corresponding symplectic leaf of M. = e*(3) with 
the cotangent bundle T*S of the unit two-dimensional Poisson sphere [7]. 

If we put C\ — (7,7) = 1, then, dividing characteristic polynomials (l47l) on det(A — Aid) 
we get the standard definitions of elliptic coordinates u, v on the sphere and their nonholonomic 
deformations u, v: 

P (\) - a + a + 7 ' _ (x-u)(\-v) 

A — a\ A — a,2 A— as (A — ai)(A — a2 j(A — 



and 



e d (X) = g( 7 ) ( r 7 1 a a-da0 + 7l(l-da 2 ) + ^-do3)- f |()j 

\ A — ai A — 02 A — 03 

(A-u)(A-v) 



(A - ai)(A - as) (A - a 3 ) ' 
respectively. Function g(7) in coordinates u, v reads as 

g(7 ) = (l-di)(l-rfv) (5Q) 

&{1 > (1 - d ai )(l - da 2 )(l - da 3 ) 1 ' 

Equation (|4"5)l is a standard definition of the elliptic coordinates u, v on the unit sphere, whereas 
equation (|4T))) determines the nonholonomic elliptic coordinates u, v. In [5], the coordinates of 
the form (|49p are used for the integration of the sphere-sphere problem under the name quasi- 
spheroconical coordinates, see discussion in [4"B] . 

We have to point out that our aim is the calculation of the variables of separation without 
any additional assumptions. Thus, we have to calculate the conjugated momenta p u ,Pv and 
p u , p v in framework of bi-hamiltonian geometry. It is easy to prove that in our case we have 
identical Stackel matrices 

S = ( ^ ^ ] and Sd ^ 



(51) 



u —v l \ — u —v 

and identical Stackel potentials 

= uH 2 -H 1 , C/ 2 (1) = vH 2 -H 1 , 

U[ 3) =uH 2 -Hi, Ujp = v H 3 - Ht , 
where 

H 1 =n i = (M, u), H 2 =H 3 = (M, M) . 

According to 44, 45], the notion of the Stackel potentials allows us to find unknown conjugated 
momenta using the Poisson brackets only. 

For instance, the following recurrence chain of the Poisson brackets 

4>x = {u, Uf'}i, 4> 2 = {u, (f>i}i, . . . , (pi = {u,<f>i-i}i (52) 
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breaks down on the third step </>3 = 0. It means that (u,p u ) is the second order polynomial 
in momentum p u and, therefore, we can define this unknown momentum in the following way 

n _ 4>i \^ Ijlmjaj -a m ){ai - u) 

Pu — -J— — / Eijm 7T7 w w 7 M-i \?o) 

02 rr 1 2(u - v)(a,j - u)(a m - u) 

13m 

up to the canonical transformations p u — > p u + f(u)- As above, Eijk is a totally skew-symmetric 
tensor. 

Similar calculation with (v,Pv) yields to the definition of the second momentum p v . In 
nonholonomic case we can perform completely identical calculations too. The results obtained 
so far can be summarized in the following statement. 

Proposition 6 The initial coordinates x = 7, M are expressed via elliptic coordinates u, v and 

Pu j Pv 



/ (u - a,i)(y - en) . . 

7* = 1/7 h r, i^j^m,. 

y (a,j - a l ){a m - di) 



(54) 



2£ijm f Yj'~fm\Q J j &m) 



((di - u)p u - (di - v)p v ^j , 



u — V 

In terms of the nonholonomic elliptic coordinates u,v and p u ,p u the same variables look like 



V^n"^ M / U 't (V ' ai) v ^^m t (55) 
(1 - du)(l - dv) U (dj - di){d m - di) 



Mj = — - — (di - u)(l - du)p n - (di - v)(l - c(v)p 1 , 

u — v V / 



where g(7) is given by [5t 

It is simple combination of definitions of the Casimir functions Ci 2, coordinates (|48H49I) . mo- 
mentum (1531) and other momenta. □ 



4.3 The reduction of the Poisson brackets 

The two sets of variables of separation u, v,p Ul p v and u, v, p K , p„ are the Darboux variables with 
respect to the brackets {., .} (f2"Tj) and {., .}d (|2"2"|) on the symplectic leaf C\ = 1 and C2 = 0. 
Of course, we can identify these variables and get the diffeomorphism cf> : M. — > A4, which 
pulls back the nonholonomic bracket {., .}d to the standard Lie-Poisson bracket {., .} on the Lie 
algebra e*(3) at C 2 = 0. 

Proposition 7 At (7, M) = the Poisson bracket {-,-}d WM between the variables 7, M 
coincides with the Lie-Poisson bracket {., .} between the variables 



7j = Jsh) (l - d(j,j)a^ jj , j 



= 1,2,3, 



(56) 



\ n^(i-d(7.7)oi) v^iw 



' 1 dVg(7)(M,A7) 7j 



T/iis mapping identifies variables u, v wf/j the usual elliptic coordinates u, v on the sphere, which 
were used by Chaplygin fllf . 
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It is a direct sequence of the previous Proposition. □ 

So, at (7, M) = we can map the nonholonomic Poisson bracket to the standard Poisson 
bracket on the cotangent bundle of the sphere. It means that any integrable system on the 
sphere has an integrable counterpart with respect to the nonholonomic bracket and vise versa. 
The list of the known integrable systems on the sphere can be found in [H [7J [2H] ■ 

In the next section we prove that we can not identify the Euler top and the nonholonomic 
Chaplygin ball using this mapping because they have different separated relations even at 
(7,M)=0. 



4.4 Separation relations 

Substituting variables 7, M (|54H55p into the Stackel potentials (|5Tj) . we obtain a pair of separa- 
tion relations (|4"0)) for the Euler top and the Chaplygin ball. These separated equations define 
some algebraic curves and we can say that the equations of motion ([1]) are linearized on the 
symmetrized product of these curves. 

Proposition 8 In holonomic case at u> = AM the variables of separation lie on two copies of 
the hyperelliptic curve of genus one 

C (1) : 4( ai -x)(a 2 -x)( a3 -x)y 2 -(xH 2 -i/ 1 )=0, (57) 

where x = u, v and y = p u ,p v . 

In nonholonomic case at u — A^M the variables of separation lie on two copies of the 
following hyperelliptic curve of genus two 

C (3) : 4(l-dx)(a 1 -x)( a2 -x)(a 3 -x)y 2 -(x# 2 - J ff 1 )=0, (58) 

where x = u, v and y = p u , p„. 

Initial variables as functions on variables of separation are given by ([541 and ([55)) . It allows 
us to express integrals of motion Hi 2 (1191) in terms of variables of separation. Substituting the 
resulting formulae for H\ 2 into the separated relation we prove this proposition. □ 

In fact, we obtain the variables of separation and the separated equations geometrically, 
i.e. without the equations of motion, the time variable and the underlying Hamiltonian or 
conformally Hamiltonian structures. We only suppose that the foliation defined by the integrals 
Hl,2 GU) on symplectic leaves of the corresponding Poisson brackets is bi-lagrangian foliation. 

However, in order to get the solutions of the separated equations x(i) and y(t) we have 
to explicitly introduce a time variable t. Solving separated equations with respect to H\ 2 one 
gets the Hamilton functions for the Euler top 

4«(oi - u)(a 2 - u)(a 3 - u) 2 4u(oi - v){a 2 - v)(a 3 - v) 2 

H l = Pu + Pv> ( 59 ) 

u — v v — u 



and for the Chaplygin ball 

Hi = -V- -,V-i ^ ~, p . ((i0) 



4v(l — du)(oi — u)(o2 — u)(oa — u) , 



4u(l - dv)(a 1 - v)(q 2 - v)(a 3 - v) 2 

~r Pu ■ 

V — U 

By definition the variables of separation are canonical variables and, therefore, we have 



Xl - x 2 x 2 - Xl 

(61) 

STJ x 4x 2 ^P d ( Xl ) 4 Xl ^P d (x 2 ) 
{H 1 ,x 1 } d = {Hi,x 2 \d= 1 ■ 

Xl - X 2 X 2 - Xl 
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Here variables xi,2 are coordinates of separation u,v or u, v, respectively. Polynomials P(x) 
and P(x) are the polynomials of degree 4 and 5 in x variable 



P(x) = (ai - x)(a 2 - x)(a 3 -x)(xil 2 - -ffi), 

P rf (x) = (l- < fx)(ai-x)(a2-x)(a 3 -x)(x J ff 2 - J ff 1 ). 
On the other hand, according to (|2U)) , the brackets (|61l) arc equal to 

2 dxi 2 2 c?xi 2 

{i?i,xi. 2 } = and {iJi,xi j2 }d = tt- 

/i dt 

where 



[i = l, and ^ = Vg(7) = 



(1 - du)(l - dv) 



(1 — dai)(l — da2)(l — ^03) 



are Jacobi multipliers (|12II13[) . 

For the Chaplygin ball, in order to get the solutions xi^fi) of the equations of motion, we 
have to consider the Jacobi inversion problem for the equations 

„ f f dxr f dx- 

/8i-2 / [i d dt = / — =^=+ ' 



(62) 

xidxi /" x 2 dx 2 



where /3i : 2 are the constants of integration. The change of time variable (|24[) reduces these 
equations to the standard Abel- Jacobi equations [TTJ [53] . 

It is easy to prove that the right hand side in fa (jUJ) coincides with an additional Euler- 
Jacobi quadrature emerged in the Jacobi last multiplier theory. Of course, for the Chaplygin 
ball this quadrature can be obtained without any change of time variable. 



5 Generalizations of the nonholonomic Chaplygin ball at 

C 2 = 0. 

Equations of motion, Poisson brackets, Lax matrices and classical r-matrices for the Chaplygin 
ball are deformation of the same objects for the Euler top by parameter d. Moreover, at C\ = 1 
and C2 = we know how to deform the corresponding variables of separation and the separated 
relations. Because at C\ = 1 and C2 = our phase space is equivalent to the cotangent bundle 
T*§ of the unit two-dimensional Poisson sphere, we can obtain similar deformations of other 
integrable systems on cotangent bundles of the Riemannian manifolds. 

5.1 The 2x2 Lax matrices. 

In variables of separation we deal with the uniform Stackel systems (|59H60j) and, therefore, we 
can get 2x2 Lax matrices associated with the Abel- Jacobi equations (|62j) in a standard way, 
see [21 [251 [371 GUI EH] as well as the relevant references therein. 

According to [37 , 38, 41 , let us introduce the following functions on the canonical variables 
of separation and spectral parameter A 

MA) = -\{ H 2 ,e(X) } h d (X) = { H 2 , e d (X) }^ , 
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and 



/(A) 



{if 2 ,MA)} ( 



4(1 - dX) 



trA-2(u + v) \ 
1 H : — I e d {X)H 2 



1 - 2dX 



1-dX 



Here e(A) and ed(X) are given by (I48H49[) . 

For the brevity below we will use the following denotation, index x is a white space for the 
Euler top and its Hamiltonian generalizations, whereas x = d for the Chaplygin ball and the 
corresponding generalizations. 



Proposition 9 At (7, M) — the Lax matrices 



h x 

fx 



-h x 



2tiL 



satisfy to the Lax equation 
d 



(J-X 

2 



^(A),^4(A) 



dt~— 2 { Hl '^} x 
Here z' = {z, H{\ x is a time derivative up to Jacobi multiplier U0il3\) . 



(63) 



(64) 



The proof is a straightforward verification of the equations (|64|) for the given Lax matrices. □ 
As usual, substituting A = x into the determinants of the Lax matrices 



which are equal to 



detJSf(A) 
detJ%(A) 



det^ x (X) = -h 2 x (X) - e x (X)f x (X) , 

XH2 — H4 
4(ai - A) (o 2 - A) (01 - A) 



4(1 - dA)(oi - A)(a 2 - A) (01 - A) ' 



one gets separated relations (f57|) and (l58l) because e x (x) = and h x (x) = y. 

In [11] Chaplygin reduces the generic case at (7, M) 7^ to the particular case at (7, M) = 
0. By applying the inverse map to the Lax matrices ([32]) one gets the Lax matrices for the 
generic case of the nonholonomic Chaplygin ball. These matrices and the corresponding r- 
matrix algebra will be studied in a forthcoming separate publication. 

Matrices ^ffe(A) are associated with the uniform Stackel systems and, therefore, they satisfy 
to the linear r-matrix algebra (f3~4")l with the well-studied dynamical r-matrices [T¥, 28, 37, 38ll4Tj. 
In contrast with the previous 3x3 Lax matrices (|32[) it allows us to obtain some well studied 
generalizations of these 2x2 matrices in the next paragraphs. 



5.2 Chaplygin ball and separable potentials. 

We are going to demonstrate that the Chaplygin ball at (7, M) = is still integrable in the 
force fields associated with a huge family of the so-called separable potentials [H Q31 US] . 

It is well known of how to get various generalizations of the separable systems using the 
deformations of their separated equations [IS] . For instance, let us consider the following 
deformations of the separation relations (|57|) and ([58)) 

4(ai - x)(a 2 - x)(a 3 - x) y 2 - (xif 2 - H{) + V(x) = 
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or 

4(1 - dx)(a x - x)(o 2 - x)(a 3 - x) y 2 - (xff s - Hi) + V(x) = , 

where the potential V is some function on x. Usually, potential V is a linear combination of 
the trivial separable potentials V m = a m x m , where m is a positive or negative integer [3] H5] . 

In order to get the same deformations in the initial variables 7, M we can use the generating 
function [?§] 

e x {\) 

or the determinant of the corresponding deformations the Lax matrix Jz? x (A) (1551) 

0(A) 



fx^ fx + 



e s (A) 



MN 



Here 0(A) is a parametric function on spectral parameter and [£(A)]mjv is a linear combination 
of the Laurent projections of £(A) by A [T4] [37l [38] . 

For example, if V = ax 2 one gets the integrable system 

ff 2 = (M,M) + a(7,A 7 ), H^M^M) ^—(^A" 1 ^, (65) 

aia 2 a 3 

which can be identified with the Neumann system on the sphere, and its nonholonomic coun- 
terpart 



H 



2 



(M, M) + ag( 7 ) ((7, A 7 ) - d(A 7 , A 7 )) , (66) 



Hi = (M, A 3 M) - a g( 7 ) ( (7 ' A 7) - daia 2 a 3 



aia 2 a 3 

Another nonholonomic analog of the Neumann system with the polynomial in 7 potential 
has been proposed by Kozlov [26] at 

rr 2 , / _, , \ , a(ai-x)(0 2 -x)(03 -x) 

y = -ax + si + a 2 + 03 x H , (67) 

1 — ax 

see integrals of motion in f|85f> . At (7, M) = this system is separable in the Chaplygin 
coordinates [17] . 

If V = /3x 3 we obtain a forth order polynomial potential on the sphere 



,( 7,A- X 7 ) 
aia 2 a 3 

and its nonholonomic analog 



Rx = (M, AiM ) + P ^' A 7) f ( 7 , A 7 ) - trA) , (68) 



Hi = (M, A 3 Af) + /9 g( 7 ) ( (7,A ' 7) - daia 2 a 3 ) x 



aia 2 a 3 



5(7) (( 7 , A 7 ) - d(A 7 , A 7 )) - trA] . (69) 



Similarly we can get other well-known integrable systems on the sphere [H HH] > such as Braden 
and Rosochatius systems, and their nonholonomic counterparts separable in the nonholonomic 
elliptic coordinates. 



5.3 Deformations of natural Hamiltonian systems on Riemannian spaces 
of constant curvature. 

At ( 7 , M) = the Euler top is a dynamical system describing free motion on the two- 
dimensional sphere, which may be identified with a particular case of the Gaudin magnet [28] . 
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It is well-known how to describe similar ./V-dimensional integrable systems on any Riemannian 
space of constant curvature and then how to add separable potentials to these systems, see 
[131 [551 133 ISFJ HOI SI] and references within. There are many different tools to investiga- 
tions of such systems (4J [34] . Below we will use only one of them based on separation of 
variables method. 

The key ingredient of this construction is 2 x 2 Lax matrix associated with the si (2) Gaudin 
magnet [HI [55J HOI [H], which is completely defined by the rational function e(A) on a given 
Riemannian manifold. The list of all admissible functions may be found in 22 . 

Let us start with elliptic coordinates (q,p) on the sphere Sat in iV + l-dimensional Euclidean 
space M.N- In this case deformation of the free motion similar to Chaplygin ball consists of three 
steps: 

• we have to change the separation relations from 

N+l 

Y[(a j -q i )p 2 i -(q?- 1 H N + --- + q i H 2 + H 1 )=0, i = l,...,N, (70) 

3=1 

to 

(1 - d Ql ) J] - ffi)P? - {Q.i~ l H N + ■■■+ q t H 2 +H 1 )=0, (71) 

3=1 

compare with (|57H58I) : 

we have to change the definition of qj in term of cartesian coordinates in M^v from 

N+1 2 ti n { \ ^ 

k =i A ak rij=i (A-OjO 

AT+l o 



that implies J2i=i "il = 1> to 

e<i(A) = g( 7 ) 2^ T - = T1 N+i n r « ( 73 ) 

k=i X ~ ak n,=i( A - a 3') 

where g( 7 ) is defined by residue of ed(A) at infinity; 

• we have to change the time variable in order to attach some nonholonomic physical mean- 
ing to the proposed pure mathematical integrals Hi , . . . , Hn and the corresponding equa- 
tions of motion, see [T9l [21] . 



Remind that rewriting separated relations (|70l) or (TTTj) in the Stackel form (|40l) 

tn 

$ « = Pi + S v{^i) H j = °> i = 1, . . . ,m, 

3=1 

we can easily determine integrals of motion in a standard way |37[ 138] 

m 

H k =Y J C k iPh 

i=i 

where C = S^ 1 is the inverse matrix to the Stackel matrix S, which is a standard transpose 
Brill-Noether matrix with entries divided by n^i 1 ( a i ~ 1i) or (1 ~ ^I^YlfJi 1 ( a j ~ Qi)- The 
same integrals may be rewritten in the following form |41j 

#*=resU =00 X N - k e-\X). 
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Of course, instead of one parametric deformation we can consider multi-parameter deformations 
replacing terms (1 — dqj) and (1 — dcij) on J|(l — d m qj) and J|(l — d m a,k) in (|7ip and (|75|) . 
respectively. 

On the other hand instead of (|72[) we can start with any other coordinate system and the 
corresponding separated equations on the Riemannian spaces of constant curvature [281 141] . 
For instance, we can take elliptic coordinates 

N , n N 



fe=l 

or parabolic coordinates 



e,; 



e(A) = A - 2 1N - ^ = -^j- (75) 

fe =i A efc IL=i (A-ei) 

in TV-dimensional Euclidean space. 

In order to consider systems with potential we can add Stackel potential U(qi) to i-th 
separated equation ([70]) and ([7T]) . 

For all these deformations we easy calculate 2x2 Lax matrices in terms of variables of 
separation, because all these deformations are uniform Stackel systems associated with various 
hyperelliptic curves [28 | l37l 138) 140 } |4T] 

So, there are not mathematical problems in the construction of such "nonholonomic" dy- 
namical systems associated with any orthogonal coordinate system and their potential gener- 
alizations. The main problem is definition of a suitable change of time variable, which may be 
associated with an interesting physical model. 



6 Generalizations of the nonholonomic Chaplygin ball at 

( 7 ,M)^0. 

In the previous section we consider various deformations of our dynamical systems at (7, M) = 
using the variables of separation method for the Hamiltonian systems and the suitable time 
reparametrization. We proceed by discussing some possible deformations of the Chaplygin ball 
in generic case. 

As above index a; is a white space for generalizations of the Euler top and x = d for 
generalizations of the Chaplygin ball. 

Let us consider the deformation of the equations ([T]) 

M = M x lo + 7 x b , 7 = 7 x cj , (76) 

where uj = A. X M and vector b is an arbitrary function on 7 and M. We want to discuss only an 
existence of integrals of motion in involution with respect to the Poisson bracket {., which is 
a necessary condition for the Liouville integrability of the corresponding Hamiltonian systems. 
The invariant measures and the corresponding time transformations are considered in review 

m- 

It is clear that Hi ,2 © remain the integrals of equations ([TBI and upon the same basis we 
can identify these integrals with Casimir functions (fT8|) on our Poisson manifold. It allows us 
to look for two additional integrals H3 and Hi in involution with respect to the same Poisson 
brackets {.,.} flST) and {.,.} d (22} 

{n 3l Hi} = 0, {H 3 ,H4} d = 0, 

where {., .}d is the "nonholonomic" deformation of canonical bracket {., .} on e*(3). For the first 
Poisson bracket all possible integrable deformations are well known J2j[7j[34]. So, we can try 
to get "nonholonomic" deformations of the Lagrange and Kowalevski tops, or of the Kirchhoff, 
Clebsh and Steklov-Lyapunov systems. 
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If the Hamilton function reads as 

2H = Ha = (M, u) + 2V(j) , lj = A X M, 
then in holonomic case equations ([75]) are identified with the Euler-Poisson equations [3J [7J 

dU dH dH f 

M = MX— +7 X— , T^7X S , H=-n i , (77) 

whereas in nonholonomic case first equation has to be replaced to 

M = MX— +7 x-, (78) 

according to the procedure of elimination of the undetermined Lagrange multipliers, see review 
[5] and references within. 

6.1 Linear integrals of motion. 

Let us briefly consider the Lagrange top [3JI7JGH] and its nonholonomic twin [IUJ[20]. Recall that 
Lagrange top is a special case of rotation of a rigid body around a fixed point in a homogeneous 
gravitational field, characterized by the following conditions: the rigid body is rotationally 
symmetric, i.e. two of its three principal moments of inertia coincide, and the fixed point lies 
on the axis of rotational symmetry. In much the same way second system is the Chaplygin ball 
with the rotationally symmetric mass distribution in the homogeneous gravitational field. 

Proposition 10 If u) = A X M, 

ai = a,2, and b =(0,0, 63), 

then the integrals of the equations |77| )- f7iSp 

H 3 = (M,M) + 2af 1 (6,7) and H4, = (M, oj) + 2(6, 7), (79) 

are in the involution with respect to the Poisson brackets {., .} \21\) and {., ,}d \22\l . respectively. 

The proof is a straightforward calculation of the Poisson brackets between integrals of motion. 
In holonomic case the linear in momenta integral 

K=(b,M)=M 3 , {H kl K} = 0, k = 1,2,3, 

can be obtained from the quadratic integrals (|79j) in a standard way 



y/y.4 - aiV.3 = \/a 3 - aiM 3 = y/a 3 - a\K. (80) 
In nonholonomic case the linear integral looks like 

^f 13 ^-^ ) ' ^ K ^ = ° k = l,2,3. (81) 



K 



It can be represented via quadratic integrals (|79p and the Casimir functions according to the 
relation 

Ui - ai H 3 + (ai - 03) (1 - dai(7, 7)) K 2 - (7, A/) 2 = . 

l-rfai(7, 7 ) 

In both cases the equations of motion (|76[) are equivalent to the Lax equation 
^ = [L,n + AB], L = M+^, A 6 E. 

(tt A 
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It is obvious, that the Lax matrix L satisfies the linear r-matrix brackets (|34|) with the same 
r- matrices (j3"5j) and ([3"rJ]l . 

In holonomic case by uj = AiM we have b x (M — lu) — at a\ = = 1 and b\ = 62 = 0. 
It allows us to get another well-known Lax representation for the Lagrange top [2j [34] : 

— = [L,0 + AB], L = AB + M + y • (82) 

In nonholonomic case b x (M — uj) 7^ and we have no such Lax matrix at all. Of course, it is 
a superficial argument because the main point is that the nonholonomic system is related with 
the genus three algebraic curve instead of the elliptic curve for the Lagrange top. 

Namely, using the Euler angles and their conjugated momenta, for the Lagrange top we 
can easily prove that the pair of canonical variables 



73 = cos(6>), Pu = — — 2 = — : Ja\i W = 1 i 

11 + 72 sm (0) 



lies on the elliptic curve defined by equation 

2 2ub + /3 — a 3 a 2 aAau + l) 2 
u z — 1 (w J — 1)^ 

where we fix the values of the integrals of motion 

Hi = irf, 1 ) = l, H 3 = (j,M) = t, K = a, H 4 = f3. (83) 
For the nonholonomic system canonical variables 

U = U, p u = —^ {u,p u } d = l, 

\/g 

satisfy to the following separated equation 

-2 , 2ub + f3 ~ a 3 p a 2 ai{au + gl) 2 

u z — 1 ^1 + dai{u 2 — l)){u 2 — l) 2 

Here we fix the values of the integrals of motion as in (|83| and 

(1 -d ai ) 2 



1 — <iai + d{a\ — a%)u 2 ' 1 + da\(u 2 — 1) 

l + rfai(u 2 -l) 

— — vs- 



If £ = (7, M) = 0, one gets the elliptic curve, but in generic case rewriting the equation (|84|) in 
polynomial form we obtain the algebraic curve of genus three. 

Only at (7, M) = we can get the solutions in the terms of elliptic functions and, therefore, 
only in this particular case we can try to reconstruct the Lax matrix associated with the elliptic 
curve. 

Three different bi-Hamiltonian structures for the Lagrange top have been obtained in 
[43] . These structures are related with different variables of separation and, therefore, different 
quadratures. If we get similar dynamical Poisson bivectors for its nonholonomic counterpart, 
one gets various quadratures, which could be associated with the distinct Lax matrices and 
underlying r-matrix algebras. 

6.2 Second order integrals of motion 

For the Kirchhoff problem, the integrable cases by Kirchhoff, Clebsch, and Steklov-Lyapunov 
are known. In this section we begin with the Clebsch case 
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Proposition 11 If u = A X M then the integrals of the equations (|77[)-J7£|) 

H 3 = (M, M) - (A 7 , 7) and U A = (M, w) + (A v 7 , 7) (85) 
are in the involution with respect to the Poisson brackets {., .} i21\) and {., ,}d (2ty)< respectively. 

As above proof is a straightforward calculation of the Poisson brackets between integrals of 
motion. 

In holonomic case we have the well-studied Clebsch problem. In nonholonomic case this 
deformation of the Chaplygin system has been proposed by Kozlov |26j in framework of the 
Euler-Jacobi integration procedure, i.e. without notion of the Poisson bracket. 

There are some different Lax matrices for the Clebsch model [3 [34]. For example, 

L = AA + M + 1^ . 

A 

We can not directly generalize this matrix to the nonholonomic case, because A^ ^ (JSJ in 
contrast with A = above. We suppose that the nonholonomic Kozlov system is related with 
the algebraic curve of higher genus and, therefore, the corresponding Lax matrices will be more 
complicated deformations of the known Lax matrices for the Clebsch problem. 

At (7, M) = the Clebsch system becomes the so-called Neumann system on the sphere, 
which is separable in the elliptic coordinates u, v (|48l) [7]. Its nonholonomic counterpart is the 
separable system in Chaplygin coordinates u, v (|4T))) [TT] and, therefore, we can get 2x2 Lax 
matrices Jz?(A) for this nonholonomic system as well. 

The Clebsch case is equivalent to the Brun case of integrability in the Euler-Poisson equa- 
tions [7] and, moreover, it is trajectory isomorphic to the Kowalevski gyrostat [23 . We can 
hope to get a nonholonomic analog of the Kowalevski top by using similar isomorphism. 

Now let us briefly discuss the integrable Steklov-Lyapunov case in the KirchhofF equation 
and the corresponding integrals of motion 

H 3 = (M, M) - 2(M, A 7 ) + (7, C 2 7 ) , 

(86) 

Ui = (M, u) + 2 (A/, A v 7 ) + (A 7 , C 2 7 ) , 

where C = diag(a2 — a.3, as — a%, a\ — CI2). These integrals are in the involution with respect to 
the first bracket {^3,^4} = 0. 

If we replace w — AM on lu = AdM in T-La then {7^3,7^4}^ = only if two parameters 
a» coincide with each other. So, for the nonholonomic bracket {., .}d we have to propose some 
more complicated deformations of the integrals of motion (1861) . 

It is known that the Steklov-Lyapunov system is equivalent to the integrable system on 
the sphere with forth order potential (|68|) |40j . We suppose that a similar transformation of 
the system (|69j) separable in nonholonomic elliptic coordinates allows us to get a nonholonomic 
counterpart of the Steklov-Lyapunov system. 

7 Conclusion 

We consider two very similar dynamical systems, which evolve on coadjoint orbits of Lie algebra 
e(3) and their non-trivial symplectic deformations. 

Close ties between the integrable Euler top and the nonholonomic Chaplygin ball allow 
us to get Lax matrices, r-matrices and bi-hamiltonian structure for this nonholonomic system. 
Moreover, in framework of the Jacobi method of separation of variables we describe a huge 
family of separable potentials, which can be added to nonholonomic Hamiltonian and briefly 
discuss how to get the iV-dimensional nonholonomic systems on the Riemannian spaces of 
constant curvature. 

In [TT] Chaplygin transforms the generic case of the rolling ball to the particular case of 
horizontal angular momentum (7, M) = 0. It allows us to solve the equations of motion using 
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the same variables of separation u, v (|49l) . which will be the non-canonical variables with respect 
to initial Poisson bracket {., .}d (l2"2")l after this map. We will discuss this Chaplygin map in 
framework of the Poisson geometry in separate publication, as well as the corresponding 2x2 
Lax matrices and the underlying r-matrix algebra. 
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